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I .  INTRODUCTION 

The  MICOM  Quiet  Radar  program  [1]  is  a  multi-year, 
exploratory  development  effort  to  build  and  test  a  short- 
range  air-defense  system  radar  with  Anti-Radiation  Missile 
(ARM)  immunity.  By  transmitting  a  low-power,  bi-phase  mod¬ 
ulated,  continuous-wave  waveform  in  conjunction  with  an 
ultra-low  sidelobes  antenna  and  a  f requency-ag *le  carrier 
frequency,  it  is  possible  to  reduce  ARM  lockon  capabilit ies 
to  ineffective  ranges. 

Previous  effort  has  determined  the  probability  of  false 
alarms  and  detections  for  the  Quiet  Radar  Processor  by 
using  Monte  Carlo  simulations.  [2]  Further  effort  was 
directed  toward  determining  the  effect  of  Constant  False 
Alarm  Rate  (CFAR)  techniques  on  processor  performance.  [3] 

The  objective  of  this  effort  was  to  determine  proba¬ 
bility  of  detection  for  a  given  false  alarm  rate  for  the 
Quiet  Radar  Processor  by  performing  covariance  matrix 
transformations.  The  analysis  included  range  cell  averag¬ 
ing  CFAR. 

Section  II  contains  the  analytical  development  used  in 
determining  processor  performance.  Section  III  presents 
the  matrices  that  describes  the  processor  elements  and  man¬ 
ipulations  required  for  a  covariance  analysis  of  the  pro¬ 
cessor.  Section  IV  presents  the  Quiet  Radar  system 
parameters  used  in  the  performance  analysis.  Section  V 
presents  performance  results  obtained  from  analysis. 

II.  PERFORMANCE  ANALYSIS 

The  block  diagram  for  the  2-D  CFAR  processor  is  shown 
in  Figure  1.  This  is  a  linear  system  up  to  the  point  where 
Z  is  calculated.  The  system  can  be  analyzed  by  a  procedure 
contained  in  a  Raytheon  report.  [4]  A  succinct  presenta¬ 
tion  of  the  analysis  follows.  The  input  X  is  represented 
as  a  column  matrix  of  the  complex  (i.e.,  I  and  Q  channels) 
input  sample  values.  It  follows  that: 


a  =  A  X  y  =  C  a  G  =  Wy  F  =  DG  .  (1-4) 
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Figure  1.  2-D  CFAR  processor  structure. 

Due  to  the  linearity,  the  output  at  F  can  be  described 
by  a  Gaussian  distribution  when  Gaussian  signals  are  the 
input  to  the  system.  Also,  the  input  can  be- separated  into 
a  sum  of  components,  viz.,  ground  clutter,  noise,  and 
target  signal.  Each  component  can  be  analyzed  separately 
using  superposition.  The  prime  objective  is  to  determine 
the  variance  at  the  FFT  output.  This  is  a  mathematically 
tractable  problem  for  the  Gaussian  ^igjials.  Let  M 
represent  the_covariance  matrix  of  X,  Mi_the  covariance 
matrix  of  a ,  M2  the  covariance  matrix  ofY,  and  M3  the 
covariance  matrix  of  G.  It  follows  that: 


-  -  -  T 

A  M  A 


(5) 


M2  =  C  M1  CT 

m3  =  w  m2  WT  (7> 


A  similar  transformation  could  be  used  to  find  the 
covariance  matrix  of  F.  However,  this  i]5  not  necessary 
because  the  variance  of  each  element  of  F  is  all  that  is 
required.  Consequently,  the  analysis  uses  the  L-point  FFT 
algorithm,  superposition  of  the  I  and  Q  signals,  and 
separation  of  the  real  and  imaginary  parts  of  the  Ffc 
element  of  F  to  obtain 


Rk 


L 

=  £ 

i,j=l 


m 


ij 


cik  cjk 


(8) 


a 


2 


Ik 


L 

£ 

i,  j=l 


m 


ij  ik  jk 


(9) 


6 


-#'"W  V-'WII. 


where  m^j  terras  are  the  elements  of  M3, 
cjk  -  ^[r-  “  1)k] 

djk  =  sin[rL  (=>  -  1)k] 


(10) 

(11) 


and  the  k  subscript  on  the  variance  represents  the  kfch 
frequency  cell  of  the  FFT  output.  Combining  the  1  and  Q 
channel  results  yields  the  real  part  of  F,  ,  i.e., 

R(Ffc)  ,  and  the  imaginary  part  of  F^,  i.e.,  I 
each  be  normal  with  variance  a  2  +  a  2,  i.e. 


<Fk> 


to 


R(Fk)  and  I  (F^)  e  N(0, 


Rk 


+  o 


Ik 


(12)  I 


This  result  holds  for  the  jth  range  bin  and  the  kfch 
frequency  cell  for  either  ground  clutter  or  noise.  A 
change  in  notation  is  used  to  represent  this  feature,  viz., 
for  noise  N 


°Njk  °NRk  +  °NIk  * 


Similar  results  hold  for  ground  clutter,  g.  Thus, 


2  2  ,  2 
a . ,  -  a„ . ,  +  a  . , 

jk  N]k  gjk 


(13) 


(14) 


The  magnitude  unit  of  Figure  1  will  change  the  Gaussian 
distribution  of  Fjk  into  an  exponential  distribution  at 
2  jk  r  i.e.. 


P  (Z 


-Z  ../2c  “ 

3 k  Dk 


(15) 
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Calculation  of  the  probability  of  false  alarm  for  a 
fixed  threshold  in  the  kth  frequency  cell  yields 


PFA 


(16) 


where 


Ybk  -  W2ojk 


(17) 


When  CFAR  techniques  are  used,  the  threshold  is  not  fixed 
but  is  a  random  variable.  It  is  possible  to  calculate  the 
expected  value  (i.e.,  average  value)  of  the  PFA  as 


PFA 


oo 

f  m  p<Ybk)dV  • 

o 


(18) 


The  density  functions  for  the  threshold  are  dependent  on 
the  CFAR  techniques. 

A  2-D  CFAR  which  averages  the  kfch  frequency  cell  of 
an  N-range-bin  window  will  have 


where  t ^  relates  the  range  bin  of  interest  to  the  range 
bins  used  in  the  CFAR  window  and  K2  is  a  threshold  con¬ 
stant  used  to  specify  a  false  alarm  probability. 
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The  development  Cor  the  probability  of  detection  fol¬ 
lows  a  similar  procedure,  i.e.. 


PDk  /  PDk  p(Ybk)dYbk 

o 


(20) 


For  a  Swerling  I  target  the  results  are 


P 


Dk 


-Ybk/(1  +  *> 


(21) 


Dk 


TkK2 


M 


for  2-D  CFAR 


+  N  ( 1  +  x) 


(22) 


where  x  is  the  s ignal-to- interference  ratio  for  the  range 
bin  and  frequency  cell  of  interest. 

At  a  given  range,  the  ground  clutter  backscatter  coef¬ 
ficients  are  assumed  to  be  constant  over  the  CFAR  window. 
However,  a  Wiebull  distribution  p(a°)  is  assumed  for  the 
range  dependency  on  these  coefficients.  The  performance 
dependency  on  a°  is  represented  as  PDk(°°) •  It 
follows  that  the  expected  value  can  be  obtained  from 


<PDk>  =  J  PDk(a°)  p(a°)  da0  .  (23) 

o 

These  procedures  were  implemented  by  Raytheon  in  a  com¬ 
puter  proqram.  The  program  is  a  hybrid  of  equation  ori¬ 
ented  calculations  and  Monte  Carlo  simulation.  The  Wiebull 
statistics  of  Equation  (22)  are  evaluated  by  Monte  Carlo 
procedures.  The  selection  of  a  2-D  CFAR  threshold,  i.e., 
threshold  from  range  bins  below  or  above  bin  of  interest, 
is  not  determined  by  Monte  Carlo  methods.  Instead,  the 
average  values  of  the  thresholds  are  determined  and  the 
largest  average  value  is  used  to  select  the  technique.  The 
mathematical  description  of  this  is  given  below.  The  range 
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bins  below  the  bin  of  interest  yield  a  2-D  CPAR  threshold 
of 

K,  -N 

\  E  Zik  (2 

K  K  2No  j=-l  DK 
-lk 

where  2-ik  represents  the  "below"  (-1)  bins  and  the 
kth  frequency  cell.  The  average  value  is 


Ak  - 


TkK2 


K2°-lk 


Similar  results  hold  for  the  2-D  CPAR  threshold  determined 
fromthe  range  bins  above  the  bin  of  interest,  i.e.,  Bk 
and  Bk.  The  program  determines 


Ybk  =  Max(V  • 


The  selection  process  is  actually  accomplished  as 


2ook2  Ybk  "  Max(K2o_lk2/  K2°lk2 


Once  the  threshold  is  selected,  then  the  results  of 
Equation  (21)  are  used  to  calculate  the  probability  of 
detection. 


COVARIANCE  TRANSFORMATION  MATRICES 


A  digital  processor  which  is  a  candidate  for  the  Quiet 
Radar  has  been  designated  as  D-8  [2]  .  A  block  diagram 
model  used  in  the  mathematical  analysis  is  given  in 
Figure  2  and  corresponding  input/output  relationships  are 
given  in  Table  1. 

The  D-8  processor  will  be  analyzed  by  the  covariance 
matrix  performance  analysis  presented  in  the  previous 
section.  It  is  readily  noticed  that  the  systems.  Figures  1 


Figure  2.  Model  used  Cor  analysis  of  configuration  D-8. 


TABLE  1.  INPUT/OUTPUT  RELATIONSHIPS  FOR  CONFIGURATION  D-8. 


MATRIX 

DIMENSION 

SYMBOL 

RELATION  TO  INPUT 

Input 

NIN  X  1 

X 

Coder 

Output 

NIN  X  1 

*C 

Bc.  X 

MTI 

Output 

M  X  1 

a 

A  •  Xc  =  ABCX 

Decoder 

Output 

M  X  1 

I 

BSc*a=  BSC*!CX 

LPF 

Output 

NFFT  X  1 

y 

C  •  8=  cbscabcx 

W  i  nd  ow 
Output 

NFFT  X  1 

cf 

W  •  T=  WCBscABcX 

Where  NIN  =  Total  No.  of  Inputs,  M  =  No.  of  MTI  Outputs, 
and  NFFT  =  No .  of  FFT  Points. 
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and  2,  are  the  same  with  the  exception  that  D-8  contains 
three  additional  elements.  They  are  (1)  coder,  (2) 
analog-to-digital  (A/D)  converter,  and  (3)  decoder. 

The  A/D  coverter  is  a  non-linear  device  and  this  type 
of  analysis  is  for  linear  systems  only.  Hence,  the  A/D 
converter  is  ignored.  For  many  cases,  the  receiver  noise 
will  dominate  over  quantization  noise.  But  if  quantization 
noise  is  sufficient  to  effect  processor  performance,  it 
could  be  applied  to  input  of  analysis  as  white  or  colored 
Gaussian  noise. 


A  discussion  on  binary  phase  coding  and  its  performance 
effect  in  the  Quiet  Radar  is  given  in  Reference  5. 
Therefore,  only  a  few  comments  concerning  the  coder  and 
decoder  are  presented.  The  coder  is  simply  a  multiplier 
used  for  coding  the  video  signal  (normally  performed  at  RF 
prior  to  transmission)  and  the  decoder  is  also  a 
multiplier  used  to  decode  the  video  signal.  The  decoder 
provides  range  cell  discrimination  when  coder  and  decoder 
have  matched  codes,  i.e.,  the  in-range  channel.  When  coder 
and  decoder  have  codes  that  are  not  matched,  the  decoder 
output  will  be  a  video  signal  modulated  by  a  shifted 
version  of  the  code,  i.e.,  the  out-of-range  channel.  Bc 
and  Bsc  are  the  coder  and  decoder  matrices  which  are  used 
in  the  covariance  analysis. 

It  can  easily  be  shown  that 


®SC  *  A  •  Bc  = 


A  for  In- Range  Channel 
A  •  T5SC  for  Out-of-Range  Channels. 


Therefore,  by 
output  covariance 


neglecting  the  out-of-range  channels, 
matrix  equation 


the 


M3  *  ( W  •  <!  •  Bsc  •A*Bc)*M*(W»<T»  Bsc  •  A  »TTC )  T 


reduces  to 


The  following  is  a  brief  description  of  the  matrices 
given  in  Table  2  for  the  covariance  analysis.  The  input 
covariance  matrix  H  is  given  as 


M 


ml,l  rol,2 
m2 , 1  m2 , 2 
•"NIN,! 


ml,NIN 

m2,NIN 

^IN  ,NIN 


This  is  a  symmetric  matrix  and  the  element  m^j  can  be 
determined  from  the  interference,  i.e. ,  noiseor  clutter 
correlation  function,  by 


m^j  =  R  (t)  t  =  (  i-j)|  .  T 


Thus,  only  NIN  elements  need  to  be  calculated. 

The  MTI  matrix  A  for  a  two-pulse  canceller  is  given  as 


A 


[lOOOO. . . 

[o][ioooo 
[oo]  [lOOOO 


[ 


-i][ 

1XNDEL 
...  -i]  [ 

...  -i][ 


i  ][ 


0 


0 

0 


10000 


] 

] 

] 


where  NDEL  is  the  number  of  input  samples  before  an  MTI 
output.  NDEL  will  be  an  integer  multiple  of  the  code 
length. 
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TABLE  2.  COVARIANCE  MATRICES  RELATIONSHIPS  FOR 
CONFIGURATION  D-8 


MATRIX 

DIMENSION 

SYMBOL 

TRANSFORMATIONS 

Input 

Covariance 

NIN  X  NIN 

M.  =  (B  *  A  *  B  )  •  M  • 

Coder 

NIN  X  NIN 

eL 

1  sc  c 

(B  •  A  •  B  )T 

sc  c 

MTI 

M  X  NIN 

A 

Decoder 

NIN  X  NIN 

Bsc 

LPF 

NFFT  X  M 

C 

M3  =  (W  *  C)  •  M1  * 

Window 

Funct ion 

NFFT  X  NFFT 

W 

(W  •  C)T 

Window 

Output 

Covariance 

NFFT  X  NFFT 

M3 

Where  NIN  =  Total  No.  of  Inputs,  M  =  No.  of  MTI  Outputs 
and  NFFT  =  No.  of  FFT  Points. 


The  LPF  matrix  C  is  given  as 


C 


[Cl  c2*  •  -CnfiLT  ]  [  0  ] 

[°J  [C1  c2  •••CnfiLt][  0  ] 

1XNSNEW 

[o]  [o]  [Cl  C2  ...CNpiLT]  [  0  ] 

[°]  .  [o][ci  C2...  Cnfilt  ] 


where  the  finite  impulse  response  filter  coefficients  are 
Ci,  C2'***CNFILT  an<3  each  row  initially  contains 
(N-l)  *  NSNEW  zeros  where  N  is  the  row  number  and  NSNEW  is 
the  sample  rate  reduction  factor. 

The  window  function  matrix  W  is  given  as 


wi  0 

0  W2 

0  W3 


w  = 


0 

0 


WNFFT 


where  the  diagonal  elements  are  the  window  function 
coef f icients. 

Since  the  input  covariance  matrix  M  is  an  NIN  x  NIN 
matrix  where  NIN  is  typically  4000  to  6000,  this  prohibits 
use  of  simple  matrix  multiples  to  perform  transformations. 
It  would  require  an  excessive  amount  of  memory  to  store 
covariance  matrix  M,  i.e.,  (4000)2  to  (6000)2  or  16  to 
36  million  words  of  memory.  Therefore,  it  is  necessary  and 
possjLble^to  calculate  one  row  of  W»C»A  and  _then  one__row  of 
J_W»C> A)«_M  and  finally  calculate  one  row  of  M3  =  (W»C»A)» 
M»(W»OA).  Then  continue  this  process  until  all  NFFT  rows 
of  M3  are  determined.  The  computer  analysis  takes 
advantage  of  this  phenomenon.  Appendix  A  contains  the 
program  listings  and  input  list. 

IV.  QUIET  RADAR  SYSTEM  PARAMETERS 

Two  different  configurations  of  the  Quiet  Radar  D-8 
Processor  will  be  studied.  The  appropriate  system 
parameters  for  each  processor  are  given  in  Table  3. 
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TABLE  3.  COMPARISON  OF  SYSTEM  PARAMETERS  FOR  PROCESSOR  #1 
AND  #2 


PARAMETERS 

PROCESSOR  #1 

PROCESSOR  *2 

Carrier  Frequency 

^c 

10  GHz 

10  GHz 

Code  Length  (PN  Code) 

NC 

31  Bits 

63  Bits 

Sample  Rate 

4  MHz 

5  MHz 

Samples/Code  Bit 

NSC 

2 

1 

MTl  Delay 

NDEL 

62 

126 

LPF  Length 

NFILT 

124  Taps 

166  Taps 

*LPF  Wait 

NWAIT 

124 

166 

Sample  Rate  Reduction 

NSNEW 

62 

83 

FFT  Length 

NFFT 

64  Points 

64  Points 

Look  Time 

1.023  msec 

1.1042  msec 

Weighting 

Hamming 

Hamming 

*LPF  wait  represents  number  of  input  samples  required 
before  LPF  output  is  used,  i.e.,  NWAIT  >  NFILT. 
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V.  PERFORMANCE  RESULTS 


Since  the  analysis  contains  several  variable 
parameters,  e.g.,  s ignal-to-noise  ratio,  clutter-to-noise 
ratio,  clutter  spread,  CFAR  window  width,  etc.,  all 
possible  performance  results  are  too  numerous  to  perform. 
Therefore,  the  analyses  were  performed  for  the  parameters 
that  have  been  used  in  previous  work  [2] ,  [3]  . 

The  performance  study  parameter  set  used  for  both 
processors  is  given  in  Table  4. 

The  results  obtained  are  plotted  in  Figure  3-14.  For 
both  processor  configurations  using  the  parameter  set  in 
Table  4,  the  probability  of  detection  for  several 
probabilities  of  false  alarm  in  each  of  the  frequency  cells 
are  shown. 

For  both  processors,  several  observations  can  be  made 
about  the  performance.  The  performance  degrades  as  the 
number  of  range  bins  in  the  CFAR  window  decreases.  This  is 
easily  explained  since  a  better  estimate  of  the  noise  is 
obtained  by  use  of  more  range  cells. 

The  loss  of  processor  performance  in  the  oth 
frequency  cell  is  due  to  attenuation  of  both  clutter  and 
target  at  oHz  by  the  MTI.  For  frequency  cells  close  to 
cell  zero,  the  clutter  residue  out  of  the  MTI  raises  the 
threshold,  thus  lowering  the  probability  of  detection. 

The  LPF  reduces  the  performance  in  the  upper  frequency 
cells  due  to  attenuation  beyond  the  target  velocities  of 
interest.  One  other  obvious  observation  is  that  as  the 
input  s ignal-to-no ise  ratio  level  is  reduced,  the  processor 
performance  decreases. 

It  is  readily  observed  that  processor  #2  performs 
better  than  processor  #1.  The  larger  dwell  time  for 
processor  #2  allows  for  more  input  samples  and,  therefore, 
the  low  pass  filter  response  of  processor  #2  can  be 
improved  over  processor  #1.  One  advantage  for  processor  #1 
is  the  hardware  reduction  possible  by  combining  the  decoder 
and  LPF  coefficient.  This  combination  is  made  possible 
since  the  number  of  LPF  coefficients  is  an  integer  multiple 
of  the  number  of  code  bits,  i.e.,  124/31  =  4. 


TABLE  4.  LIST  OF  PERFORMANCE  STUDY  PARAMETERS 


PARAMETER 

VALUE 

S/N  (Input  Signal-to-Noise) 

-21,  -24  dB 

C/N  (Input  Clutter-to-Noise) 

46  dB 

erf  (Clutter  Spectral  Width) 

8  Hz 

CFAR  Window  Width 

4,  8,  16  Range  Bins 

PfA 

IQ-3,  10“4,  10“5,  10-6 
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S/N  a  -21  dB,C/N  =46  dB ,  8Hz  GROUND  CLUTTER 
HAMMING  WEIGHTS,  4  RANGE  BINS  IN  CFAR 
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S/N  =  -21  dB,C/N,  =  46  dB,  8Hz  GROUND  CLUTTER 
HAMMING  WEIGHTS,  16  RANGE  BINS  IN  CPAR 
WINDOW,  2-PULSE  CANCELLER,  124  TAP  FIR 
FILTER,  64  POINT  FFT 
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S/N  =  -21  dB,C/N  »  46  dB,  8Hz  GROUND  CLUTTER 
HAMMING  WEIGHTS,  4  RANGE  BINS  IN  CFAR 
WINDOW,  2-PULSE  CANCELLER,  166  TAP  FIR 
FILTER,  64  POINT  FFT 
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S/N  =  -21  dB , C/N  =  46  dB,  8Hz  GROUND  CLUTTER 
HAMMING  WEIGHTS,  16  RANGE  BINS  IN  CFAR 
WINDOW,  2-PULSE  CANCELLER,  166  TAP  FIR 
FILTER,  64  POINT  FFT 
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Figure  12.  Detection  performance  of  a  124-tap  LPF  D-a 

processor  with  S/N  =  -24  dB  and  a  4  range  bin 
CFAR. 


S/N  »  -24  dB,C/N  =  46  dB,  8Hz  G HOUND  CLUTTER 
HAMMING  WEIGHTS,  8  RANGE  BINS  IN  CEAR 
WINDOW,  2-PULSE  CANCELLER,  124  TAP  FIR 
FILTER,  64  POINT  FFT 
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Figure  13.  Detection  performance  of  a  124-tap  LPF  D-8 

processor  with  S/N  =  -24  dB  and  an  8  range  bin 
CFAR. 


S/N  -  -24  dB,C/N  =  46  dB ,  8Hz  GROUND  CLUTTER 
HAMMING  WEIGHTS,  16  RANGE  BINS  IN  CPAR 
WINDOW,  2-PULSE  CANCELLER,  124  TAP  FIR 
FILTER,  64  POINT  FFT 


Noiioaiaa  jo  uinavaowd 


Figure  14.  Detection  performance  of  a  124-tap  LPF  D-8 

processor  with  S/N  =  -24  dB  and  a  16  range  bin 


APPEND 

PROGRAM 
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The  development  of  the  Quiet  Radar  Processor  covariance 
analysis  was  performed  on  the  AP-120B  array  processor  and  a 
PDP-11/10  host  computer  contained  in  the  Data  Acquisition 
and  Analysis  System  [6] . 

The  simulation  was  written  in  Fortran  utilizing  the 
high-speed  processing  capability  of  the  AP-120B.  The 
program  contained  calls  to  AP  math  library  subroutines  and 
to  AP  assembly  language  programs  which  simulate  processor 
elements  [2].  The  input  to  the  program  is  the  same  as  the 
Monte  Carlo  simulation  input  in  Reference  2.  Hence, 
providing  capability  of  ready  comparison  between  Monte 
Carlo  and  covariance  analysis  results. 
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FRECtUUG  fAvS  £LmmK«N0T  ¥IlM£> 


CARD  NO 


List  o£  Inputs 
VARIABLES 


FORMAT 


NC  =  Length  of  pseudo-random  code 
NSCLK  =  Number  of  samples  per  code  clock 
period 

NFILT  *  Number  of  FIR  filter- 
samples 

NSNEW  *  Sampling  rate  reduction  factor 
NHAIT  *  Delay  in  FIR  filter  output 
NP  =  Number  of  two  pulse  canceller 

stages  1615 

NDEL  *  Number  of  samples  in  MTI  delay 
NFFT  =  Number  of  FFT  samples 
NCINT  =  Number  of  non-coherent 
integration  samples 

NFA  =  Number  of  different  thresholds 


s imulated 

IC(I),  1*1,...,  NC  =  One  period  of 
code 

(consists  only  of  ones  and  zeros) 

the 

8011 

FSAMP  * 
VMAX  * 

NBIT  * 

Sampling  rate 

A/D  converter  saturation 
voltage 

Number  of  bits  in  A/D  converter 

2E15.4 

15 

NOISDB  * 
SIGDB  * 
FDOP 

KT 

Noise  power  in  dB 

Target  return  power  in  dB 

Target  doppler  frequency  in  Hz 
Target  delay  in  number  of 
samples 

3E15.3 

15 

DISCDB  * 

SIGMAF  * 
XM 

A 

FIXCDB  * 
KF 

Distributed  clutter  power 
in  dB 

Clutter  spectral  spread  in  Hz 
Clutter  DC-to-AC  power 
rat  io,m2 

WEIBULL  parameter 

Fixed  clutter  power  in  dB 

Fixed  clutter  delay  in  number 
of  samples 

5E15.3 

15 

List  of  Inputs  (Concluded) 


CARD  NO. 

VARIABLES 

FORMAT 

6-8 

Gaussian  Curve  fit  parameters 
[X(j)l,  l  y  ( j  )  ]  and  [  z  ( j  )  ]  for  j*l, 

•  •  •  r  6 

6E12.5 

9-33  * 

FIR  impulse  response  samples. 

NFILT  numbers,  five  per  card. 

5E16.8 

34-36* 

Pre-FFT  weighting  sequence.  NFFT 
numbers,  five  per  card. 

5E16.8 

47 

ALP(I)  1=1 , . . . ,NFA  =  Multiplication 
factor  to  control  the  threshold 
levels 

5E16.8 

48 

N RUN ( I ) ,  1=1,..., NFA  =  Number  of 
Monte-Carlo  trials  for  each 
threshold  setting. 

1615 

49  * 

CFARK(I) ,1=1, . . . ,NFA  =  Multipli¬ 
cation  Factors  to  Control  Range 

Cell  Averaging  CFAR  Threshold  Level 

5E16.8 

50 

CFAFK (I), 1=1,... ,NFA  =  Multipli¬ 
cation  Factors  to  Control  Frequency 
Cell  Averaging  CFAR  Threshold  Level 

5E16.8 

•Denotes  inputs  used  by  covariance  analysis. 


35 


C  FORTRAN  KKUGKAM  EUR  PERFORMANCE  ANALYSIS 

C  UK  UU1ET  KAUAR  SIGNAL  PROCESSOR  HY 

C  COVARIANCE  MATRIX  TECHNIQUES.  THE 

C  MATKIX  SlAfcS  REQUIRED  USE  UK  AN  AP-120 

C  An  D  A  PUP-11.  THE  MATKIX  TRANSFORMATION 

C  IS  UK  THE  KORM  IKwCA) ,M. (E ftCA )**T  rHERE 


c 

M  IS 

THE  INPUT  COVARIANCE  MATRIX 

N  In 

X 

NiN 

c 

A  IS 

CANCELEK  MATRIX 

M 

X 

NlN 

c 

c  IS 

LUft  PASS  E lLlEK  MATRIX 

NFF1 

X 

M 

c 

ft  IS 

EFT  WEIGHTS  MATRIX 

NETT 

X 

NFF1 

c 

E  IS 

KPT  MATRIX 

NKFT 

X 

NFFT 

c 

WHERE  n f  t'  1 

*  NU.  UK  FE1  PUInTS 

c 

M 

s  (NEE i-1 JNSnEr+NFILT  s  INPUTS  TO 

LPF 

c 

mn 

*  wSC*NC4l\  *  MfNUELsTUT,  NU.  UE  INPUTS 

c 

NE  1LT 

s  NU.  OE  E1R  KILTER  CuLE  F  TC1ENTS 

c 

hSNEW 

s  sampling  hate  reduction  factor 

c 

M<eL 

s  NU.  OE  INPUTS  bEEUKL  MTi  OUTPUT 

c 

NSC 

S  Input  samples  per  code 

c 

NC 

S  CUDE  LENGTH 

c 

A 

s  NU.  UE  TRANSMISSIONS 

C  TYPICAL  NUNHEKS 

C  NUT  *  h4 

C  M  =  AO  30 

C  NlN  s  <40*2 

C  MULT  s  124 

C  NSNEr  *  b2 

C  NULL  a  b2 

C  NSC  »  2 

C  NC  s  jl 

C  P  =  bb 

C  ASSUMES  DEC'UDlNG  CUMPLtT'EU 

CuMMUn/ARkA  I/Mb4  ),CHlbOU),CEAK*(y),FTLT(2UO),SH(64), 

♦  SU.CS(b4)  ,.SIGNSl«»4)  ,SCAbt.T(b4J 
COMMUn/CUNST/NC.N>  1LT  »NSNEft*Np»NUE.L*NEEI  »NiN»NE'A»T»NWAlT  i 

♦  NUTS,SJu,OTSCLT ,nCE  E i # A.SIGMAE  ,  Twl ,N 
REAL  MUlS 

CALI  TURK 
CALL  APRUN 
Call  si Gnu 
CALI  Oak 
s  i  u  p 

t>  i  * 

SUPKUU  1  ll<t  TURK 

C'-M»m</AP*Mi/ft  (b4l  «CK  (bUO)  ,CE  ARM  (*  J  ,E  TL  j  (200)  .Si  l  (b4)  » 

*  SiLLS(b4j , Si GNS(b4) .SCALE! (b4) 

CUr  (■  un/Uihb  l  /  i«C  ,  NE  I  L  l  ,  NSNE*  ,  NP  ,  NULL  ,  NP  K  1  ,  N  IN  .  nK  a  ,  T  ,  NwA  I  T  t 

*  hOTS.STG.UISCLT  t  NCKK  T  »  A  .SIGN AE  ,lNT,N 
UT.'EwSTUh  *ETGHl(b4) 
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JHIS  PAG®  IS  BEST  QUALITY  PRACTICABLE 
Jggyi  OOI Y  FULiiiiStiEiJ  TO  ODC 


ulMKNSltm  Ail'll  1  0  U  J  *  A.  t<  U IV  l  V  )  .CfrArMO)  ,  1  <_  (  1  u  o  J 

Ht'AL  AiUio ,  >«u ihCih 

t^uiVALtMt.  I  A  IN  ,  Abl’  »  iM<On  #  Cr  Air  ,  1C)  ,  l  *  (  *  M  on  1  ) 

1  ij  HO  1 

ktALi(b,l)  '4C  » i»dCLit  «  lit-  iLi  *  N£n  f. »  ,  i*  n  A  1  ]  ,  t  , 

1  hin*hUM#Ni» 

KtA0(b,2  )  l  iCIU  *1S1  MC) 
i  t UkaAI 1 1 bl b  j 
2  KOhAiA'I  lb  01  1  ) 

J  KUHwAl I2tlb. 4,  lb) 

KfcAb  l  b  #  J  J  r  bA*l  ,  V  MAA  ,  Nt*ll 
Kf.Aulb,  4  Jnui&iie  ,oiiiL>b  ,t  Out',  M 
4  1 dkn Al  I3blb..l. lb) 

t\LMl>(b/bJ  i.’ii'Cbt'  #i>ib  hmi-  #ai,a,H  IaUM,k  t 
b  I*  tlx*  A 1 1  bt  l  b ,  i  $  1  b ) 

x  t AO ( b  #  b  J  l AlN(l) ,1=13# ioj 
b  iuxMAi Ibti 2,b  1 

KbAL.  (b.  7  J  lML.ll  1  J  ,  i  =  l  ,1.1-  I LJ  ) 

7  LJxrflllbtlo.bJ 

HtAblb./)  l*tllinUU,ls),ni'H) 

HtAMb#  7  MmOM  1 ) ,  1  =  1  ,nt  A) 
K6;AL>lb#»HMKUNlu,lai,M.A) 
b  J-OWwAillblb) 

KtAolb,  7  HO  AKM1)  ,1*1  ,iM  A) 

KtAblb.l)  n 

K  t  A  U  ( b  #  7HCMFMl),i*t,NM) 

KbAMb.U  n1 
L)U  12  UI.nih 
SiOCbl  I  )=0, 

61(.i«6l  11=11. 

12  CurTiMifc 

uu  17  i=l,t>oO 
O(l)=0. 

1/  C'JMlNOfc 

Uu  10  l  =  l,'«Mlil 
CKll  )=ML1  ill 
10  Corn  Ihlit 

1  =  1  ,/lSAHF 

»•  AL  1=0.2  J02bb 

u U.= !  .  /I  i>AMP 

NO  1  1=MJ*N1  te  l 

NOlb  =  k.AH(lACl*^UifaUU) 

biG=fc.AMlHCl*SfOL<H) 

UiSCL'l  =lc.Al'  1 1  AC’l  *l)lSCua  ) 

■M  1  N 1  =N1  M 

NlN=(NlNl«l)ANSbtl»  +  i'iwAlI+Nl>*NLItL 
NHUL)  =  AU0(N1N,I»C) 

IKlNMUH.Nt.O)  Mnsn1n  +  mC»M»0L) 

iiCnLib=Mw/oC 

UAbl.L  =  OtL.*KLUAl  IMA) 
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* HIS  PAiib  IS  BjSSI  QUAim  P8ACI1CABL* 

rtwM  OOPX  i'UiwsISHJsa  ro  DDC  - 


OUTPUT 

PRINT  40,LRtLL 

FORMAT  1/ 10* , 'PROCESSOR  U-6*/ 

1  IUX.'OrELL  11ME  *  ',tl2.bl 
FRInT  11  ,NC,NSCLK«NP.NLEL,NFtT  ,N,tSAMP 
tUkMATl/lOX  ,  'COLE  PtkluL  a  ',15/ 

1  10X,'NU.  Ut  SAMPLES  PER  CLOCK  PEHlUL  =  ',15/ 


l  1uX,'nU.  of  PULSES  CAnCeLLEL  In  mTI  a  ',1b/ 
l  1UX*'  NO  ■  Ut  SAMPLES  IN  MIX  DELAY  *  ',15/ 

1  lux,  ' nu ,  Of  ffl  SAMPLES  *  ' *  lb/ 

'  lUX,'wU«  OF  RANGE  CELLS  IN  (.'FAN  RJNUUR  a 
>  IUX , 'SAMPLING  RATE  AT  INPUT  a  '  ,Ei2.b,'  HZ 
PRIM  1  J,NSNE*,NRAIf  ,NFXLX 

FUR*  Al  (  /  10X  *  'SAMPLING  NATE  RELUCTIUN  fACTUN  a 


'  ,1b/ 
HZ'  ) 


1  1UX,  '  NUMBER  OF  IKANSIEN'I  SAMPLES  DELETED 

2  IUX,  'NO,  Of  FIR  1MPULES  RfcSPUNSt  SAMPLES 
PRInT  2J 

HJNMA'l  1/10X,  'FILTER  COEFFICIENTS'/ ) 

PN1NI  24,  (t  ILIU),1*1,NF1L1  ) 

tURMAT 1 bx,btlb.b) 

PRINT  14,SlGLb,NUiSUb 
tUNWATl/lUX,  'TARGET  RETURN  plJREk  a  ',E12,b,' 
1  1UX,'NUISE  Purer  a  ' , t 1 2 .b , '  Lb'l 
PKInI  lb  ,01 SCOb ,S1GMAF , A 

FORMA 1(/10A,  'D1STK1HUTE0  CLUTTER  POrER  a  », 

1  IUX, ' CLUT 1 ER  SPECTRAL  SPREAL  *  ',E12,b,' 

2  IUX, 'RtibULL  PARAMETER  a  *,E12.b) 

PRINT  26 

FORMAK/lUX,  'rETGHTInG  CUEt  F  ICIENTS  '  /  ) 

PRIM  lb,  C  NEIGH  1  (I  )  ,Isl,NFFT) 
FURMATlbX,bElb,b) 
ap  clear 
CALL  APCLR 

CALL  VCLR(U, 1,32/67) 

CALL  APrK 
RETURN 
END 


=  ',1b/ 


, El 2  .  b  ,  ' 
HE'/ 


SUbHUUIINE  APRUN 

COM RUN /ARRAY/ A (64  ) ,CR(600),CFAKK(V),FlLT{2UU),STl(b4), 

♦  SlGCS(b4l ,SIGNS(64 ) ,SCALET164) 

C0MM0N/CUN5I/NC,«tlLl  ,NSNE»,NP,NDtL,Nt  FT,MN,Nt>  A,  1,  Nr  AIT, 

*  nuis,sig,uisclt,ncffi ,A,SIGMAF , INI , N 

DIMENSION  M(bOUO) 

HEAL  M , NUTS 

UATA  M/141, ,bV99*U,/ 


THIS  PAG£  (.  :  TJAL1M  PRACTICABLE 

5*KuM  0-ux  X  PuxCvlDXiXiiC  10  Oi->C  —  "«*•"  ~~ 


nnr: 


C  AF  INITIALIZATION 

Fi*3, 141592L54 

GKNLHATlNU  uNt  KUW  ui  CA  mackia 
i»F  .  Lfc. ,  2 

No.  UK  CUL&  •  tOUAL  ivF*NLKL*  M1  I  i<  A 
NDKL2*N0fcL*  2 
N0KL1=NLLL+1 
Uli  6  0  J*1#nF 
LmJ  bO  i  =  ivLtL  1  ,  NUKL2 
CF 1 1  +  NL>fcL2  )  s*Ch  ( 1  ♦ivUfe.L  ) 

CK(ifiHik.L)a-CHCl)+CKU4wUfc.l,) 

Ck (1 >  =C  h ( 1 J -CK t i-NUFLi 
b  0  C  u  •»  1  1 1«  0 1 
bO  CUn]  l  -vut 
C  A  F  b  A  1  A  tt-  Ini 

'')C>'  =  fJUtL*w^  +  i^  1  Li 
•vCUK  =  /OUU+  «  J  iM-i 

.vC LI  A  ^SM  Uh  +  'ii»"l 

1 1*  2  =  2  *  H  1  «.  -  1 

C  LOuF  t  Oh  GbNFKAlllvL  CLU I  IKK  AUi>  AHkaXS 

ULI  1  (/  0  1 IV  1st  ,2 

C  Of  I«  t.K  AT  J  iv  c.,  Klin  Ut-  lrtt  M  i*  Al  K  i  X 

if  I  liVl  .Fw.  I  )  Gi.  1U  10 
AFL=/.  »r  i*olLHAi  *1 
uii  0  1=1  , iv  I  n 

•*  (  l  J  =  l  l-(Ahb*U-lX  )**2)/2.) 
o  Cu«  i  |  >v  o  f. 

Ju  Cut  J  i l U r. 

C  A  L  Ij  A  F  C  L  h 

Call  vcl*  t  o  ,  i  »  32 /t>  /  j 
Call  a  f  * k 

Cali.  A  F  f  U 1 1 IV ,  i  , 1 1  1  ,  i  J 
CmlL  AFFOilCK#  l  0  o  #  M.’ v  » 2  ) 

CALL  AFFU’J  l  r-  ,  uCLK  ,  ;«  i  f-  ,  i  ) 

call  ak»l 

CALL  VMuv  l  .-CCA  l  ,*1 , 2UUU  ,  1  ,  J  1  >  ) 

I  r  1 1'vl  .to  .  /  J  CALl  vtXK/imu,  1  ,  /  o  wu  ,  |  , .  i  >  f  j 
C  A  L  L  A  F  *  A 
CALL  AKt  A  I. 

Call  a K r  fr  i 
CALL  Ai  *»h 
ioo  CUhlluOt 
F t 1  o  K  iv 
hivL 

vSOLAUlJ'l  INI;  AFFAL 

CUM" UN  /  Akh  M  I  /  i  l  t>H  )  ,Cf  (  MlO  )  ,(  I-  Af  ►  IV  )  ,  t  i  L  i  l  /UU  t  ,  <>  1  l  (  h  4  J  , 

♦  S  i  GC»>  l  *3  4  )  ,  l>  1 L.  I'oloi  )  iMAU  I  lf.fi 
CuFkUU/CCM»i>  1  /  Hi  ,  Ivf  l  L  I  *  iv.iiJ*- .  *  ,  A  l  , ..  H  L  ,  !«r  \r  I  ,  t.  i  L  »  ivf  A  ,  i  ,  I»  A  A  i  l  , 


39 


,,SBBSiaUALlW  thAUii° 
*his  PAGE  de>c  - - 

jaB»00Wl»*lSJiwlU 


\  w(llS,SlG,OiSCLl  #NCfc  f  T  #  A#SlGAAfc 
WbAL  NU1S 
f«CMsNUKL*NP-»NFlL‘i 
AP  bXbCUlb 

LUUP  2U  GbNbtfAlbb  M  J*(  MCA)  »M.  OCA)**l 
Hi  IS  AN  M-fcl  X  f*H  l  NAlklX 

Mi  IS  S'lUKbU  HUm  Ot  KUw  fc  HUM  SIAKllNL  AUOKtSS  240OU 

ltPs24UOU 

141=2000 

UU  20  1*  =  1,«FH 

LOUP  JO  GbNbKA'i bS  Uwb  RU*  OF  M  *he.Kb  M1=1»CA),m 
(JNfc.  K(jm  HAS  N1N  tLbMbNlS 
UU  30  1MU=1#N1N 

lbs  1 400U  + »Ml-l 

Caul  VNutl  luo,  1 ,  1*,  l , suu,  i  ,  NC*») 

CAbl-  S  v  b  l  ouu  ,  i  ,  1  i  00  ,  nCm  ) 

Call  vS*uli 1* , 1 , 1100, Ifc , I , 1 ) 
call  #t'«K 
iu  CUN  1  IN  llfc, 

illSlCHNOIt* 

1  b=  1  4C'00 

ljUC’fc  to  ufc  u  r.H  a  t  t.S  UNfc  Hum  t.ifc  SJ  MibKb  m i=  1  *Ca  )  , M  ,  ( *CA )  * ♦! 
U.*fc_  All*  MM. -5  Nfcfcl  t  L  b  i*!  fc.  N  1  s 
lu  iv  i * 2= r# i.fc  fc  t 

UMLl  (M)L(  lOO  ,  1  ,  Ifc  ,  1  fbOU*  1  »  N  t  H) 

Call  svt. i cul ,  i ,  l  luo , nci* ) 

Call  vslllii**^#i#  i  ioo»lt,.fc,»l#  1) 

C  A  lj  o  l4‘  *  h 

L  r,S  i  t  ♦  HOI»fc  * 

1  r,fc  =  l  e  fc  fci«fc  fc  i 

lu  OlL'jiMJf 

It  fc'  =  2  4  U  l<  L  ♦  i  « 

2 b  C  ii«  1  i  ini', 

K  K  i  I'M' 
b  ■  <  1 1 

SMbHUU  I'lufc,  Mbfcfcl 

C‘  l«"HL  r>  /  AH  K  A  I  /*  t  t:  4  )  ,CK  l  KUO  )  ,  Cfc  AKK  l  H  )  ,  fc  1  b'i  (  200  J  ,  SI  1  1 0  4  )  , 

+  S  ILLS  l  t<*t  )  ,  SlLiio  l  to-*  J  t  SC  A  Lfc  i  (  04  ) 

Ci.l"  >i  it.  /  I  i,ii.»>  i  /  iv  C  #•*1’  1  L  1  t  r»si«b*  t  Nb  »  NUfc.L  •  NF  F 1  t  n  1 W  $  NfA  »  1  f  iK  A  A  1  T  t 

|  ■* » i 1  o  0  S  I  L  f  U 1 S  C  l*  l  #  •»  b  fc  fc  I  >  A  #  »SI  Cj^i  At*  §  1  N  1  f  N 

Ml  A  I  *  *  »-> 

Hr. Mfc  ui>  «  l  i«b  r  ,MJ,  fc 

L  .j  v-  ►  lit.  t  l.  1  r.Kfc  ►  Kt  NCfc  A  fc  1  fc,  H  fcfcl 

l.t  I  S  =  i>  fc  fc  i 

I.t,'  rthll'A.  fc  MA  I  Hi  K 
r  IS  u  I.fc  fc  l  A  I.fc  fc  1  MATk  1  * 

40 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

IBDM  COPY  FUKiilStibO  TO  DsDC  _ 


n  o  r;  n 


C  F  IS  SFPAkAlbU  inlO  i>ik  A  HO  CCS  i-*AikJC t.f* 

NCUS«0 

nSlksbOOG 

ZbkU*t). 

CALL  APPUl(Ztk<J.JiibQ»I«4J 

Call  a^*u 
Fib  *l*.*Pl 
uO  uu  isi(M>lS 
0*i»i 

FinC*Hn  ♦ j 

CALL  APPU'J'IF  ikC  $  S4l*>  i  1 1  $  /  ) 

CALL  Aino 

Call  vkamf  c  j/itot  in  /b  1 1  iioco,  i  ,*L'i  o> 

CALL  VCOMmuuu#l#NCuS*MFM  ,«HS) 

CALL  WMwl  WOOL,  I  #<v.»iM#NF>l  ,  Li  i.S  ) 

CaLL  AFwk 
NCUSsNCOStl 
«<SlN*nSihF  1 
oU  CONliNub 

GbnbkAl ino  t«*i 
PbkFUKkinG  F.AJ 

CALL  kMUL ( U  #  1 » 44UUU  * 1 » i ZULU  *  I  ,!•►>!  ,  nf  ►  j  ,  .-.r  F 1 1 
CALL  PAULlbUOU,  t ,*4  0uy , I , lkULu, J ,NFF I r X,nFFi  ) 
CALL  MT*AkStU,l,0#l,ntfri#Htfr  i) 

CALL  NTKAWS(4>UUU»i(bUOU»l#wFF  i(HFF  i  ) 

CALL  APkfN 

ObibNMlniNL  dIAlomaL  bbtftbt  IS  lit 
SiUkb  tLbKfcNlS  At  AUOMFSS  Z  4  0  0  U 
MC*0 
*S*bOUO 
iO*Z4UUO 
‘4FTC*12UOO 
NFTSsiHOOM 
LO  1  OU 

CALL  VMULlnF  iC#l»HS,Mt,i  ,*C#  t  ,NPlSi 
Call  vaaim- is#kPis**s,i,MC#l*MS»i»kPiM 

CALL  SVb(A*,  1,  lUtkPIS) 

CALL  APwh 
MC*ACtNFlS 
MSsnStNFTS 
Iu*lu*l 
NF1C*NF1C-»1 
NFlS*nFiS4 1 
1UU  CUMikUb 

iFUM.bU.U  CALL  APGbi  ISiGkJj  ,  ZaOUO  ,  wl  F  T  ,  Z  J 
1 F  ( I  k  i  « bu »d i  CALL  APGbi tbIGCS « Z aOOO , *t  F  I *  Z  ) 

CALL  APnL 

HblUPh 

t«0 


41 


Iti****1'*** 

DDDC 


o  O  b  K L'  1 1 1'<  t,  C 1  M  I*  , 

Cu*u<,lju /ak  h  m  i  /  *  itu  )  ,  Cut  t>uU  )  ,C  1' Ak'M  9  )  ,11  Li  1  ^uu )  i  £  ( t><* ) , 

*  *  1  w(  ,S(f>4J  ,£lC.Ni>(04),i3CALll(b47 

Cu*  M."'.  /CO’«b  1  /  riC  »  M  1  L  i  #  Nktitb*  .  "1  ,  MliFL  ,  ht  F  1  ,  ft  ill  ,  Ml  A  ,  1 ,  *«»  A  i.  1  , 

*  M‘iO  , is  lb  r  1*1  SCL-  1  »  »  C 1 1  1  ,  A  ,b  luwAr  ,  1 M  ,  f< 

Uli'tt.SlON  *t  I  (  ILLb  J  ,  1  1  11 t>4)  ,  frl  i  I(b4) 

Ul*  t  *J>  J  b*  I-  u 1  tj  &  J  ,  hLbi’U  l  b4  ) 

Kt*.  A  t-  i-Uli* 

CALL  A*;>1CM3, 'LMlCI- A*.»  bl  '  ,0,  *»«t*  1  J 

Ue.1  Iwl.  1  i  i.r  i  ( i  by  ,  i /»  ,  u  #  u  J  i 
CAl  L  *Mibl  t  uli»CL  1  r  A  ,  *CL  ) 

I  C  l  =  1 

I I  *  S  1  I J  u 

9  LUw  (  1  M-t 

OU  it!  lC=l,NfM 
AMr'  =  Ct  maa  1  1C  7  *t* 

K 1  A  a  1 

DU  bub  1  K  =  I  ,  'if  1  1 
.''OM  HMsu.u 
but*  I  Cui-  i  1  '<Ot 

UO  bl<  jU  N  *  S  1  ,  14  ft  * 

DU  bu>4w  r-i  ,  wl r  t 

III  UJS‘U.U>  ft  7*iilbCbif»  J1WU  lo*i>lbMMK  J 
IrilsilllM 

iHI  K>  I  .  Hw  , 0 . 7  1  tM *1 .  Ul  - 1 0 

Sill  lK7a(aib*l»CAbfcIlF  7  )  /  1  CiST 

FUl*  >sl  ./(  l.KAMV  l*M  t.*£>ril(F  ))  7  7*** 

Jj'J M I  UHO*i>U«l>p(*  7+FblK  7. 

b»J40  cuillwbt 
bOju  CUiUli.UK 


DU  UlbU  IS  1,  Mr  11 
i»U  «H>  1 1  )  SiHi^kU  1  i  J  /”>  ** 
olbu  CUwlJNUb 


1  OD 

1U1 


1U41 

IDi 

lu 


7 


ft  K  1 1 1. 1 1>  *  1 U  0  7 

1  t.lK  A  A  1  1  !M,b4A,  •  LU 1 bl  K  AT)  AK  L-»  Kl.ftUblO*  7 
**1  11  l  O,  Id  7  K|AM|11A 

1  lift  k  M  1  (//IDA,  *  NO  »  Ul  Cbt.LJj  IN  Cl  AK  wlNLUft  a 
1  1uX,'C1mK  1  Hnfc  bCAIJMi  CUNbiAM  =  *  , blk.b/ 

t  1UA  ,  '  KhUBA*' 1  Li  1  1  ul  lALisl  AbAKM  s  *,bl2,57 
*  K 1 1 1  Ibilm  J 

I'UKkA  I  (  /  X  /  /  b4  A  »  '  1'HUDMbi  LI  11  Ul  l>lTbC'ii.UN,//7 

*K1  1  ir.lb  ,  10  i  )  U  .bU-Kuil  >  ,  1*1  ,M  r  I  ) 

luiNAllb4A,lb,bA,t.W.b7 

«Kl  lbi i' 111  7  IbDhlDl  11 )  ,  11*1  ,MM  7 

1 C 1* IC1 ♦ 1 


Cell'll  1  ■«  Ob 

nc.au  l  b ,  /  )1C1AHM1J#1*1  ,nFA7 
Kb AD l  b  ,  1  7* 

lFlM.Nb.U)  CO  lu  9 

FUKKAI Ibk  lb.b  J 


42 


Ob »- 

1 


»AGi* 


>  -si 


^XICASI* 


1 


tufcMA'I  ( lblb) 

KLTUHN 
fcND 

SUbKOUTlNfc  NbULL l DISC bl *  A  #  *CL ) 

OiNtNSlUN  «CL(1000) 

DU  10  1*1*1000 
K*K AN  10*0) 

»*CL(l)*(ALUG(l./ll.-K)>>»*A 
«CL(l)*01SCLl»*CLli)/2. 

10  CUNliNUb 
SUA*0. 

UU  20  1*1  *  1000 

SU#*5U**nCL(1  ) 

20  CONllNUt. 

5UM*MJM/  lUOU, 

OU  iO  1*1*1000 
WCL(1  )*01SCL1 **CL ( 1  )/60M 
iO  CUMlNUt 

NKlTfc.  (  6  *  bO  )  SUM 

bO  t  0  K  M  A I 1/bX '  wtAN  OF  UNSCALtO  CLOlltK  ChOSS  SfcCilON  StU  lb  't20. 

HtlUHN 
tNU 


c 


c 


c 


SiibKUUl 1 Nfc  S1GNU 

COMMON /ANNA  i /«(0*J,CklbOOJ,CrAKMO).f  11*1  (  200  ),Slllb4), 

♦  SlGCS(b4) *SlGNb(b41 *SCALb1 (b4J 
CuMMuN/CONSl/NC.Nl 111  * NSNfc* »  NP  *  M/fc 0 « w«- f 1  * N 1 « * N* A  * 1  * w» A IT  * 

♦  NulS*SIG*OlSCLl#NCt»l*A*SlUMAt  ,1M,N 
KbAL  nUIS 

OiHtNblUN  11(04) 

NnA11*NNA11'*1 
P1*J, 141b92bb 

CUMPUlt  FWfcwUfcNC»  lNCNl-.mt.Nlb 

r  ACi  *1  •  /  (  nkk  ltftbNb*) 


2tNu*o, 

1n11 lALHAl ION 

CALL  APPUlUtHO,  J27bto,l,2) 

call  appokhli, Jo/5o,NriLi *2) 

CALL  APPUl  («*,  JubOO,NrFT,2) 

CAuL  APftU 
NU*32UOO 
AH*2.*Pl*t  ACI 

LOOP  KUH  SIGNAL  PuwfcN  CALCULA1 1UN 
00  lb  N*l»NtPT 
AKG* AN  * ( A • 1  ) 

CALL  APPU'l(AMG»J2)N)*l*2) 

CALL  AP«0 

CALL  VNAMPC 32?00* 32 /to)* 12000  *  1 *NlN) 


43 


c 


lb 


10 

lul 

1UU 


1  Bb 


CALL  VSiN(12000,l#0,l,*iL) 

CALL  VCUS(12000,l,6000,l»MN) 

CALL  *11(0, NOEL, 6P#KIN)  , 

CALL  *11(6000, 6LEL#NP«NlL) 

CALL  F1KT( 0,30 750, 20400, NF1LT , NSNEw ,N*AIT ,NF* T ) 
CALL  FIKK600U, 30750, 26401, *KiLI,wS*t»,W*AlI,«f FT) 
CALL  MtiGhK  20400, 30600, Nfr  FT) 

CALL  CFFT6( 20400, 29000, MFFT,i) 

CALL  C  V*  ACS  (  29000, 2, 29600  il,  NtF'i  ) 

CALL  MAXV(2960Q,1,NU,NFF1) 

CALL  AP*K 
WiJsNU+1 


CUM1MJE 

CALL  APGtl (SCALE! , 12000 ,  NFFT  ,  2 
CALL  AFWL 

SIGNAL, NUlSE  ALP  CLUlfEK  AKHA1 
LJ  10  1*1,  »FU 

KA'i  i*AbS  (  S 1  bCS  (11  J/AtsS  (Si  GCS  (  1 
IF  IK All, LI . 1 ,E-J)  S1GCS (11*0, 
CUMlNLfc 
*Kllt(6,101) 

FUKAAT ( 1H 1 / // 1  OX  » 6HSI6N Ab , JK A , 
«  k 11 1 ( b  » 100)  (SCALE! (1 ) » SIGN SI 
KIK*Al(13A,tlO.H,2/A,hlb.O,2/X 


) 


UUlPUlS 

)) 


bnMJlSE,i7A,7HCLUTTt.K/) 

1),S1GCS(1),1*1,NFF1) 

,E16.0J 


Du  1  Kb  N*1,*FF1 

11 (A JsuiSCLl*SlGC6  (lO+NL  AS*  SIGNS  (E  ) 


ItS 1*1 1  IK ) 

IF (1  ESI ,(0,0,1  ltSl*i.E-10 
61 AU  )*(S1G*SCALE‘1  (A)  1/1ES1 
CUMlwut 


Kt'l'OF  i'i 
Fwu 


.  vNCVS  Xs 
.  ••• 


C  oe- 


. ,  '.A1  • 


44 


APPENDIX  B 

RANGE  CELL  AVERAGING  CFAR  THRESHOLD 


The  range  cell  averaging  CFAR  threshold  is  calculated 


N 

VTk  K  N  3?1  zjk 


(B-l) 


where  V-j-jt  is  estimated  threshold 

Zj](  is  FFT  square  law  output 
N  is  CFAR  window  length 

K  is  scale  factor  for  a  specified  average  probability 
of  false  alarm 
j  is  range  bin  index 
k  is  frequency  cell  index. 


For  a  system  as  shown  in  Figure  B-l,  if  I  and  Q  are 
Gaussian,  then  Zjk  is  an  exponential "distribution.  Thus, 


P(Z 


-*jl/2oJk 


(B-2) 


Figure  B-l.  Detection  System. 

It  can  be  shown  that,  by  knowing  Equation  (2) ,  a 
probability  density  function  can  be  derived  for  the 
threshold,  obtained  in  Equation  (1) ,  and  an  average  PFA 
can  be  found ,  i . e . , 


where 

Tk  = 


variance  of  range  cells  of  estimate 
variance  of  range  cell  of  interest 
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1 


Assume  all  range  cells  have  identical  noise,  i,e, 
or  "jk2  *  0ok2'  then 


-N 


=  N(  /fPS  -  1) 


For  example,  let  PFA  *  10“6  and  N  ■  4,  then 
K  =  4  (-41/T^6_1)  «  122.5 


CFAR  threshold  constants  used  in  the  analysis  are  given 
in  Table  B-l. 
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TABLE  B-I.  CFAK  THRESHOLD  CONSTANTS 


>FA 

N 

K/N 

K 

r3 

4 

4.62 

18 .49 

8 

1.37 

10.97 

16 

0.54 

8.64 

,-4 

4 

9.00 

36.00 

8 

2.162 

17.30 

16 

0.778 

12.45 

-5 

4 

16.78 

67.13 

8 

3.22 

25.74 

16 

1.05 

16.86 

-6 

4 

30.62 

122.05 

8 

4.62 

36.99 

16 

1.37 

21.94 
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